Zhang et at. Plant Methods 2014, 10:18 
http://www.plantnnethods.conn/content/10/1/1 



PLANT METHODS 



METHODOLOGY Open Access 



Protocol: a beginner's guide to the analysis of 
RNA-directed DNA methylation in plants 

Huiming Zhang^" Kai Tang\ Bangshing Wang\ Cheng-Guo Duan\ Zhaobo Lang^ and Jian-Kang Zhu^'^"" 



Abstract 

Background: DNA methylation is a conserved epigenetic marl< tliat controls genome stability, development and 
environmental responses in many eukaryotes. DNA methylation can be guided by non-coding RNAs that include small 
interfering RNAs and scaffold RNAs. Although measurement of DNA methylation and regulatory non-coding RNAs is 
desirable for many biologists who are interested in exploring epigenetic regulation in their areas, conventional 
methods have limitations and are technically challenging. For instance, traditional siRNA detection through RNA 
hybridization requires relatively large amount of small RNAs and involves radioactive isotopes. An alternative approach 
is RT-qPCR that employs stem loop primers during reverse transcription; however, it requires a prerequisite that the 
exact sequences of siRNAs should be known. 

Results: By using the model organism Arabidopsis thaliana, we developed an easy-to-follow, integrative procedure for 
time-efficient, quantitative measurement of DNA methylation, small interfering RNAs, and scaffold RNAs. Starting with 
simplified nucleic acid manipulation, we examined DNA methylation levels by using Chop PCR (methylation-sensitive 
enzyme digestion followed by PCR), which allowed for fast screening for DNA methylation mutants without the need 
of transgenic reporters. We deployed a simple bioinformatics method for mining published small RNA databases, in 
order to obtain the nucleotide (nt) sequences of individual 24nt siRNAs within the regions of interest. The protocol of 
commercial TaqMan Small RNA Assay was subsequently optimized for reliable quantitative detection of individual 
siRNAs. We used nested qPCR to quantif/ scaffold RNAs that are of low abundance and without Poly-A tails. In addition, 
nuclei fraction enables separation of chromatin-associated scaffold RNAs from their cognate non-scaffold transcripts 
that have been released from chromatin. 

Conclusions: We have developed a procedure for quantitative investigations on nucleic acids that are core 
components of RNA-directed DNA methylation. Our results not only demonstrated the efficacy of this procedure, but 
also provide lists of methylation-sensitive restriction enzymes, novel DNA methylation marker loci, and related siRNA 
sequences, all of which can be valuable for future epigenetic studies. Importantly, step-by-step protocols are provided 
in details such that the approaches can be easily followed by biologists with little experience in epigenetics. 

Keywords: Epigenetics, RNA-directed DNA methylation, siRNA, Scaffold RNA, Pol IV, Pol V, Nuclei fractionation 



Background 

DNA methylation is a major epigenetic mark that can 
confer transcriptional silencing of genes and transposable 
elements. RNA-directed DNA methylation (RdDM), which 
involves complimentary pairing between small RNAs 
[including small interfering RNAs (siRNAs) and piwi- 
associated RNAs (piRNAs)] and scaffold RNAs, has been 
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established in plants and mammals [1-5]. DNA methylation 
is subject to dynamic regulation during establishment, 
maintenance, and removal in response to developmental 
and environmental cues [6]. In order to understand epigen- 
etic regulation, it is important to characterize the histone 
modification patterns and to determine the levels of DNA 
methylation and small RNAs, as well as of scaffold RNAs 
that are chromatin-associated long non-coding RNAs 
(Figure lA). This study focuses on assays for DNA methyla- 
tion, small RNAs and scaffold RNAs, since methods for 
characterizing histone modification patterns such as chro- 
matin immunoprecipitation (ChIP) have been covered 
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Figure 1 An integrative procedure for studying nucleic acids in RdDM. [A] A model of RdDM in Arabidopsis. Briefly, Pol IV initiates siRNA 
production by producing single-stranded non-coding RNAs, while Pol V synthesizes scaffold RNAs that recruits Argonaute-associated siRNAs. The 
complimentary pairing between the two types of non-coding RNAs, as well as physical interactions between proteins in the methylation complex, 
results in targeting of the methyltransferase DRM2 to RdDM loci. For simplicity, not all known RdDM components are shown. [B] Chop PCR-based 
screen for mutants showing abnormal DNA methylation patterns. Red asterisks represent enzyme restriction sites. Cytosine methylation overlapping 
restriction sites will protect DNA from methylation-sensitive digestion (MSD), while enzymatic cleavage results in failure of subsequent PCR amplification 
of the target sequences. Note that the enzyme McrBc is special in that it cleaves methylated- but not unmethylated DNA sequences (see also 
Additional file 2: Table SI). [C] RT-qPCR detection of siRNAs by using stem-loop primers during reverse transcription. [D] Schematic procedure of 
nuclei fractionation, which separates chromatin-associated RNAs from RNAs that are not associated with chromatin. See Methods for details. 



extensively [7,8]. Genome-wide DNA methylation and 
siRNA accumulation patterns can be precisely profiled 
by next-generation-sequencing. At individual loci, conven- 
tional methods for obtaining data on DNA methylation and 
non-coding RNAs include bisulfite sequencing for DNA 
methylation, RNA hybridization for siRNAs, and RT-PCR 
(reverse transcription followed by polymerase chain reac- 
tion) for scaffold RNAs [9,10]. Although commonly used, 
these methods have limitations. Bisulfite sequencing is 
costly and time-consuming, and the assay requires consid- 
erable experience in several important steps of the experi- 
mental procedure including designing suitable primers. 
Northern blotting of siRNA requires large amounts of small 
RNA input and involves radioactive isotopes. RT-PCR is 



semi-quantitative and, if applied alone, does not distinguish 
scaffold RNAs from those transcripts that have been re- 
leased from chromatin. To overcome these limitations, we 
used the model organism Arabidopsis thaliana to develop 
PCR- or qPCR-based methods for studying nucleic acids in 
RdDM. These optimized approaches form an integrative 
procedure that can be followed by biologists with little ex- 
perience in epigenetic studies (Figure 1). 

Results and discussion 

The first part of this procedure is to detect alteration in 
DNA methylation levels. Alternative to bisulfite sequen- 
cing, DNA methylation levels in different samples can be 
compared by Chop PCR, in which methylation-sensitive 
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enzyme digestion of DNA is followed by PGR amplifica- 
tion of the target sequence (Figure IB). By using Chop 
PGR, we developed a method for epigenetic mutant 
screen that does not need a transgenic reporter gene. A 
large population of Arabidopsis mutants is commercially 
available from seed stock sources such as ABRC 
(Arabidopsis Biological Resource Center) and can be 
used for Chop PCR-based mutant screen. Genomic 
DNA samples were prepared by a simple home-made 
recipe that is suitable for simultaneous DNA extraction 
from multiple samples. We used the methylation-sensitive 
restriction enzyme Hae III and examined DNA methyla- 
tion at AtSNl, a known RdDM target locus. Hundreds of 
Arabidopsis mutants were screened, resulting in discovery 
of several DNA hypomethylation mutants including 
atrrp6ll alleles that were recently characterized [11]. 
In addition to Hae III, we have examined some other 
methylation-sensitive restriction enzymes and correspond- 
ingly demonstrated the efficacy of these enzymes in 
examining several RdDM marker loci in Arabidopsis 
(Additional file 1: Figure SI; Additional file 2: Table SI). 
Different restriction enzymes can be used to assay 
methyaltion in the different sequence contexts, i.e., CG, 
CHG, and CHH (H represents A, T, or C). Although Chop 
PCR examines less cytosine than bisulfite sequencing, 
observation of DNA hypomethylation by Chop PCR in 
multiple marker loci is sufficient to demonstrate an 
epigenetic function. Quantification of the difference in 
DNA methylation levels can be achieved by Chop qPCR 
(Additional file 1: Figure SIB). Parallel amplification of the 
marker locus from non-digested DNA samples is 
commonly used to show loading controls. However, we 
recommend amplifying DNA sequences without re- 
striction sites, from the same digested DNA samples, 
as internal controls, because experimental accuracy 
can be increased while the amount of required DNA 
can be reduced. Compared to traditional forward genetic 
screens that involve transgenic reporter genes. Chop 
PCR-based screen does not require generation of a 
transgenic reporter line and the subsequent mutated 
population, neither required is the process of map-based 
cloning that can be time-consuming. 

The second part of this procedure is quantification 
of siRNAs (Figure IC). The levels of siRNAs can be 
semi-quantitatively detected by RNA hybridization, 
which typically uses >20 (ig RNAs that are enriched 
in sizes smaller than 200 bp [12]. The whole process 
takes several days from RNA extraction to final results. 
Alternative to RNA hybridization, small RNAs can be 
detected by RT-qPCR that employs stem loop primers 
during reverse transcription [13]. This approach takes 
only a few hours starting from RNA extraction and 
involves the patented TaqMan Small RNA Assays. 
However, the protocol of TaqMan Small RNA Assays, 



which was initially implemented to detect miRNAs, 
does not readily enable detection of lower abundance 
small RNAs such as 24 nt siRNAs that mediate RdDM. In 
addition, these PCR-based quantitative measurements 
require, as a prerequisite, knowing the exact sequences of 
small RNAs; meanwhile mining published small RNA 
datasets for such information typically requires skills 
in bioinformatics. Therefore we developed a simple 
bioinformatics method for selecting siRNAs to quantify 
and optimized the TaqMan Small RNA Assays protocol 
for detecting 24 nt siRNAs. 

To obtain nucleotide sequences for 24 nt siRNAs from 
the regions of interest, we searched small RNA sequencing 
databases. This can be done by BLAST searches of 
published small RNA databases, such as Zhang et al. [14], 
for siRNAs that match the regions of interest. However, 
such analysis requires skills in bioinformatics. We thus 
provide an alternative analysis approach that can be 
applied in any laboratory, by using Arabidopsis AtSNl locus 
as an example. AtSNl locates within 1 kb downstream of 
the gene At3g44006 [9]. Browsing At3g44006 and its 
downstream region in IGB (Integrated Genome Brower) 
identified AtSNl as At3TE63860, which is a transposable 
element that starts at 15794808 and ends at 15794619. 
Arabidopsis small RNA nucleotide sequence information 
with genomic coordinates is available in the UDelaware 
Small RNA database [15]. Since 24 nt siRNAs are proposed 
to pair with scaffold RNAs [1-4], we screened A small 
RNAs for 24 nt siRNAs that are annotated as Waston 
strand, i.e., template DNA strand. This means that these 
siRNAs have the same sequences as the template 
DNA strand (despite the U/T difference) and thus 
can potentially pair with scaffold RNAs. To increase 
the chance of detecting individual siRNAs, these 24 nt 
siRNAs were subjected to a second screen for those 
that showed abundance of > 1 hit per milUon reads, 
according to the UDelaware Small RNA database, in 
the wild type Col-0 ecotype. Among these candidates, 
one siRNA exhibited sequence homology to a 43 nt 
probe that had been used for detecting AtSNl siRNAs 
by Northern blots [10]. Thus we chose this siRNA for 
customer TaqMan Small RNA Assay design. 

The manufacture-recommended protocol of TaqMan 
Small RNA Assay suggests using 1-10 ng total RNA 
per reaction. We were able to observe stress-induced 
regulation of miRNAs by using > 50 ng total RNA as 
input (Additional file 1: Figure S2). However, most 
24 nt siRNAs that we examined were undetectable 
under similar condition. We found that 400 ng small 
RNA (the <200 bp fraction of total RNA), instead of 
total RNA, per RT reaction is sufficient for consistently 
detecting all 24 siRNAs that we tested, including the 
AtSNl siRNA (Figure 2A). We used snoRlOl, a small 
nucleolar RNA, as the internal control to ensure equal 
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Figure 2 Quantitative measurement of 24nt siRNAs and scaffold RNAs in Arabidopsis. [A] Quantification of individual 24nt siRNAs by 
TaqMan small RNA assays. The levels of siRNAs were normalized by using snoRlOl as internal control. [B] Quantification of Pol V-dependent 
scaffold RNAs by gene/scaffold RNA-specific RT-qPCR. Actin2 was used as internal control. Error bars indicate SD, n = 3 technical replicates. All experiments 
were performed with at least two biological repeats. Representative images are shown. 



RNA input During the RT reactions, we mixed RNA with 
2 (iL 20X (manufacture-defined concentration) siRNA RT 
primer together with 2 \iL snoRlOl RT primer, instead of 
using 3 (iL siRNA RT primer alone as suggested by 
the manual. Subsequently, the levels of siRNA and 
snoRlOl were independently measured by quantitative 
PGR. Because siRNA and snoRlOl are simultaneously 
reverse transcribed in the same tube, the level of 
snoRlOl reflects RT efficiency of each sample, allowing 
for elimination of false results that may be caused by 
different RT efficiency between samples. 

The effectiveness of using TaqMan Assays to quantify 
24nt siRNA levels was confirmed by using the Arabidopsis 
mutants nrpdl-3 that is defective in RNA polymerase Pol 
IV, which is required for biogenesis of almost all 24nt 
siRNAs [14], and nrpel-11 that is defective in Pol V, 
of which mutation decreases the levels of some Pol 
IV-dependent siRNAs [16]. Compared to wild type plants 
(Col-0), nrpdl-3 exhibited depletion of all examined 
siRNAs including AtSNl, siRlOOS, solo LTR, 5TE27040, 
IGN5B, and Cluster 2 (Figure 2A). The nrpel-11 mutant 
also showed substantial reduction in the levels of the 
examined siRNAs except Cluster 2 (Figure 2A), which 
is known to be Pol IV-dependent and Pol V-independent 
[17]. The patterns observed in our RT-qPCR measure- 
ments are consistent with published RNA hybridization 
results [10,17] and whole genome small RNA sequen- 
cing results [18]. Sequences of these individual 24nt 
siRNAs are listed in Additional file 2: Table S2. Designing 
siRNA probes requires only submitting individual siRNA 
sequences to the manufacturer website of the patented 
TaqMan Assays. 

The third part of studying RdDM nucleic acids is to 
quantitatively measure chromatin-associated scaffold RNAs. 
In Arabidopsis, Pol V transcription generates scaffold RNAs 



that are of low abundance and are devoid of poly- A tail [9]. 
Therefore we used scaffold RNA-specific primers during 
reverse transcription for the first-strand cDNA synthesis, 
and subsequently performed nested qPCR. Such an 
approach successfully measures Pol V-transcribed scaffold 
RNAs, as demonstrated by comparison between Col-0 
and nrpel-11 (Figure 2B). Although the levels of scaffold 
RNAs are in proportion to the pool of total Pol V 
transcripts, sometimes nuclei fractionation is needed to 
distinguish chromatin-associated scaffold RNAs from 
those Pol V transcripts that have been released from 
chromatin (Figure ID), as was the case in our recent 
study of AtRRP6Ll that functions in retention of Pol V 
transcripts at chromatin [11]. Nuclei can be isolated with 
Honda buffer and subsequently lysed in the presence 
of urea and detergent, followed by centrifugation that 
sediments chromatin-associated RNAs in the pellet 
(chromatin-associated fraction), while leaving released 
RNAs in the supernatant (chromatin-free fraction) 
(Figure ID). As an indication of successful nuclei 
fractionation, genomic DNA should be nearly absent in 
the chromatin-free fractions (Additional file 1: Figure S3). 

To provide an easy-to-follow procedure, experimental 
details are described step by step in the Methods. 

Conclusions 

To summarize, we have developed an integrative pro- 
cedure for studying nucleic acids that are central 
players in RNA-directed DNA methylation and other 
epigenetic regulation (Figure 1). The results not only 
demonstrated the efficacy of this procedure (Figure 2; 
Additional file 1: Figure S1-S3), but also provide a 
list of novel RdDM marker loci (Additional file 2: 
Table SI) and the related methylation-sensitive restriction 
enzymes (Additional file 2: Table SI), as well as the 
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nucleotide sequences of a list of individual 24nt siRNA 
(Additional file 2: Table S2), which can be valuable for 
future epigenetic studies. Moreover, this procedure requires 
little previous experimental experience in epigenetics or 
bioinformatics, and thus can broadly benefit the scientific 
community. 

Methods 

I. Chop-PCR/qPCR 

<l-1> Preparation of genomic DNA samples 

1. Prepare 1.5 mL eppendorf tubes, each containing 3 
Chrome steel beads (BioSpec Products, Inc, 
Cat.# 11079132C). 

* Note: Steps 1 to 3 are optimized for preparing 
multiple samples. Start from Step 4 if samples are 
homogenized by other methods. 

2. Take one or two mature leaves and put into labeled 
tubes. Make sure the tubes are tightly capped. 

3. Freeze tubes in liquid nitrogen shortly, then shake 
tubes vigorously in a box to shatter the samples. 

*Note: Be aware that a capped tube with liquid 
nitrogen in it will pop open or even explode under 
room temperature. 

4. Add 500 (iL 2% CTAB (Cetyltrimethyl ammonium 
bromide) solution to each tube (or per -100 mg 
sample if homogenized by other methods), then 
incubate at 65°C for 30-60 min. 

* Note: To make 100 mL CTAB solution 
(modified from [19]): 2g CTAB, 10 mL 1 M 
Tris-HCl (pH = 75), 28 mL 5 M NaCl, 4 mL 0.5 M 
EDTA (pH = 8.0), H2O. 

5. Add 500 (iL chloroform, then votex vigorously and 
centrifuge at > 12000 g at room temperature for 10 min. 

6. Transfer aqueous phase into new tubes, then mix 
with equal volume isopropanol. 

7. Store at -20°C for > 30 min. 

8. Centrifuge at > 12000 g at room temperature for 10 
min. 

9. Wash the pellet with 0.5 mL 70% ethanol. 

10. Air dry pellets and dissolve in 50 [iL nuclease-free 
water (we use autoclaved double-distilled water). 

11. Proceed to steps in < 1-2 > for Chop-PCR if doing 
mutant screening. Alternatively, for Chop-qPCR, 
continue with DNA purification as described below. 

* Note: Although DNA samples extracted by this 
way have RNA contamination, they can be used 
for rough estimations of DNA methylation levels 
during initial mutant screen. RNase treatments 
should be performed for Chop qPCR and for 
semi-quantitative Chop PCR that is used for the 
purposes other than rough estimation. 

12. Add 2 (iL RNase A (lOmg/mL) per 500 ^L DNA 
solution, mix well and incubate at 37°C for 1 hr. 



* Note 1: Concentration of DNA solutions can 
vary depending on the yield of crude extraction 
using CTAB solution. 

* Note 2: If making RNase A solution from 
powders, boiling the solution for 15-20 minutes 
usually eliminate most of its DNases. 

13. Add equal volume phenol: chloroform: isoamyl 
alcohol (24:1:1, pH 8.0), vortex vigorously for 5 sec, 
then centrifuge at > 12000 g at room temperature 
for 10 min. 

14. (Optional) Transfer the aqueous phase to a new 
tube, repeat Step 2. 

15. Transfer the aqueous phase to a new tube, add 
equal volume phenol: chloroform, vortex vigorously 
for 5 sec, then centrifuge at > 12000 g at room 
temperature for 8 min. 

16. Transfer the aqueous phase to a new tube, add 2 
volume ethanol, 0.1 volume NaOAc (pH 5.2, 3 M), 
mix well. 

17. Store at -20°C for > 30 min. 

18. Centrifuge at > 12000g at room temperature for 
10 min. 

19. Wash the pellet 2X with 0.5 mL 70% ethanol. 

20. Air dry pellets and dissolve in nuclease-free water. 

</-2 > Enzymatic DNA digestion and PCR/qPCR 

1. Measure genomic DNA concentration using 
Nanodrop. 

2. Digest equal amount (typically 0.5 to 1 (ig) of 
genomic DNA with a methylation-sensitive enzyme 
in a 20 [il reaction mixture. 

* Note: While following manufacture protocols to 
set up reactions for each enzyme, the digestion 
time can be extended to overnight. 

3. After digestion, check DNA methylation levels at the 
enzyme-targeted marker locus by qPCR or PCR. For 
semi-quantitative PCR,use 1 (iL of the digested DNA as 
template in a 20 (il reaction mixture that contains 0.15 
unit of commercial Taq DNA polymerase. Run a DNA 
gel (1.2% agarose) to compare the PCR results from 
different samples. For qPCR, prepare 20 (il reaction 
mixture following manufacture protocols and use 1 (il 
of the digested DNA as template per reaction. 

* Note 1: DNA input amounts in different samples 
should be normalized by amplifying an internal 
control sequence that does not contain the 
restriction site of the enzyme used. Alternatively, 
non-digested DNA can be used as PCR/qPCR 
templates for input normalization. 

* Note 2: In semi-quantitative PCR, too many PCR 
cycles will probably yield in saturated amplification 
and thereby mask small differences in template 
inputs. For optimized observation of differences. 



Zhang et at. Plant Methods 2014, 10:18 
http://www.plantnnethods.conn/content/10/1/18 



Page 6 of 9 



PGR can be repeated with increased or decreased 
cycle numbers after initial assessment. A positive 
control, e.g., the Pol IV mutant nrpdl-3, should be 
examined in parallel to the wild type. 

Optional: For Chop-qPCR, follow steps 4 to 10 to 
analyze the results 

4. Export qPCR results to Excel file. 

5. Calculate average Ct value for the non-digested 
control in each sample. 

6. Calculate "ACt = Ct of the examined amplicons - 
average Ct of the non-digested contror. 

7. Calculate average ACt value of one sample 
(e.g., wild type plants). 

8. Calculate the comparative Ct values (AACt) for 
other samples (e.g., mutant plants). AACt = ACt 
(mutant) " average ACt (^na type)- 

9. Calculate relative abundance of the amplicons by 
using this function: 2^(-AACt)". 

10. Calculate average of relative abundance of the 
amplicons; and calculate standard deviation by using 
the function "=STDEV(values to be calculated)". 



* Note: After removing the bulk isopropanol 
solution, centrifuge briefly and use pipette to 
remove the residual liquid. By doing this, RNA 
pellets dissolved without drying perform well in 
the downstream RT-PCR reactions. 

9. Dissolve the pellet in 20-40 |iL RNase-free water, 
without drying, in RNase-free water. Keep on ice and 
proceed to reverse transcription, or store at -80°C. 

<//-2> Quantitative measurement of siRNAs by RT-qPCR 

1. Search small RNA databases to determine sequences of 
siRNAs to be examined (see main text for an example). 

* Note: Additional file 2: Table S2 lists some 
Arabidopsis 24nt siRNAs that have been detected 
by RT-qPCR in this study. 

2. Submit siRNA sequences to Life Technologies for 
TaqMan Small RNA Assay design by following the 
online instructions. 

3. Place the order after the designing has been 
completed. 

* Note: we suggest downloading the TaqMan Small 
RNA Assay protocol at this point, as well as to 
order the TaqMan MicroRNA Reverse 
Transcription Kit and the TaqMan Universal PCR 
Master Mix as recommended in the protocol. 

Use TaqMan MicroRNA Reverse Transcription Kit 
to prepare reverse transcription (RT) master mixture 
and keep on ice. The following recipe is for one 
reaction: 
100 mM dNTPs 0.15 
MultiScribe TM Reverse Transcriptase 
(50 unit/^L) 1.0 

lOX Reverse Transcription Buffer 1.5 [iL 
RNaseOut (40 unit/^L) 0.9 
RNase-free Water 4.26 |iL 
Total Volume: 7 (iL 
In an RNase-free PCR tube, add 400 ng small RNA 
in a volume of 5 [iL. Then add 2 (iL snoRlOl RT 
primer and 2 (iL siRNA RT primer. Mix well gently. 
400ng small RNA is sufficient for consistently 
detecting all 24 nt siRNAs that we have examined. 
For certain abundant siRNAs such as siRlOOS that 
originates from 5S rDNA repeats, 100 ng small 
RNA is sufficient. 
Incubate the RNA-primers mixture at 85°C for 
5 min, followed by 60°C for 5 min. Then put the 
mixture on ice immediately and for at least 1 min. 
Add 7 (iL RT mixture to each RNA-primers mixture. 
Mix well gently. Keep in ice for 5 min. 
Perform RT reaction using the following conditions 
consecutively: 16°C for 30 min, 42°C for 30 min, 
85°C for 5 min. Keep the samples on ice after RT 
reaction. 



II. SiRNA RT-qPCR 
<//-7> Isolation of small RNAs 

4. 

1. Homogenize <100 mg plant tissues in liquid 
nitrogen by using mortar and pestle. Transfer the 
powder into a 1.5 mL eppendorf tube (pre-chilled in 
liquid nitrogen) and add 1 mL RNAzol® RT 
(Molecular Research Center Inc) 

2. Add 0.4 mL water per 1 mL of RNAzoP RT. Mix 
well and store at room temperature for 5-15 min. 

* Note: Autoclaved distilled water is good enough 
for this step. 

3. Centrifuge the mixture > 12,000 g at room 
temperature for 15 min. 5. 

* Note: this step precipitates DNA, protein, and 
polysaccharide. 

4. Transfer the supernatant to a new tube and add 
0.4 mL 75% ethanol per 1 mL of supernatant. Mix 
well and store at room temperature for 10 min. 

5. Centrifuge > 12,000 g for 8 min at room temperature 
or at 4°C. 

Note: this step pellets RNA > 200 bp. 6. 

6. Transfer the supernatant to a new tube and add 0.8 
volume of isopropanol. Mix well and store at 4°C for 
30 min. 7. 

* Note: the mixture can be stored at -20°C for 
longer time. 8. 

7. Centrifuge >12,000g for 15 min at room 
temperature or at 4°C. 

8. Wash the pellet with 70% isopropanol twice. 
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* Note: Because TaqMan stem-loop primer 
recognizes a specific short sequence and the 
qPCR reverse primer is complementary to the 
stem-loop RT primers (see Figure IC in main text), 
genomic DNA contamination, if any, will not be a 
problem for this experiment. Thus no-RT control is 
unnecessary. In fact, we did not detect qPCR signals 
even using genomic DNA. 

9. Prepare the qPCR mixture as the following: 

TaqMan Small RNA Assay (20x) 1.0 

cDNA from RT reaction (Step 8) 1.2 

TaqMan Universal PGR Master Mix II (2x) 10.0 |iL 

Nuclease-free water 7.8 |aL 

Total Volume: 20 |iL 

* Note 1: run qPCR individually for snoRlOl and 
the siRNA and perform 3 technical repeats. 

Note 2: The patented TaqMan Small RNA 
Assays provide qPCR primers in the tubes labeled 
as "TaqMan Small RNA Assay (20x)" without 
giving primers information. 

10. Perform qPCR in a real-time PGR instrument using 
the following program: 

Step I: 50°C (optional) 2 min 

Step II: 95°C (Enzyme activation) 10 min 

Step III: 95°C (Template denaturing) 15 sec 

Step IV: 60° C (Anneal/extend) 1 min 

Step V: go to Step III, 39 times (totally 40 cycles) 

Step VI: End 

* Note: TaqMan Small RNA Assay PGR signals 
will be detected in the FAM channel by the qPCR 
instrument. 

11. Export results to Excel file. 

12. Calculate average Ct value for snoRlOl in different 
replicates. 

13. Calculate "ACt = Ct of siRNA - average Ct of 
snoRlOr from the same RT reaction. 

14. Calculate average ACt value of one RNA sample 
(e.g., wild type plants). 

15. Calculate the comparative Ct values (AACt) for 
other samples (e.g., mutant plants). AACt = ACt 
(mutant) " average ACt (wiia type)- 

16. Calculate relative siRNA abundance by using this 
function: 2^(-AACt)". 

17. Calculate average of relative siRNA abundance; 
and calculate standard deviation by using the 
function "=STDEV(values to be calculated)". 

* Note: Mean value and standard deviation can 
also be calculated based on multiple biological 
replicates. The expression patterns, e.g., release of 
transcriptional silencing, should always be the 
same for different biological replicates, although 
the standard deviation between different biological 
replicates is usually greater than that of technical 
repeats from one single biological replicate. 



III. Scaffold RNA RT-qPCR 

<///-/> Quantitative measurement of scaffold RNAs without 
nuclei fractionation 

1. Extract total RNA by using Trizol (Invitrogen) 
following the manufacturer's protocol. 

* Note: Scaffold RNAs may be a few hundred 
nucleotides, so for this experiment we do not 
recommend using those commercial RNA 
extraction kits that preferentially extract 
RNA > 200 bp. 

2. Remove DNA contamination in the RNA samples by 
using TURBO DNA-free™ (Ambion). 

* Note: We typically treat 10 |ig RNA with 
TURBO DNase in a 20 (iL reaction, and use 
maximum amount of RNA during RT reaction 
while following the manufacturer RT protocol. 

3. Perform RT reaction by using Superscript™ III 
First-Strand Synthesis SuperMix (Invitrogen). In 
each reaction, use mixed gene/scaffold RNA-specific 
primers (the scaffold RNA and the internal control 
ACTIN2). For a 10 (iL reaction, use 0.5 (iL pre-mixed 
primers (2 mM each). 

* Note 1: We have tried using some other RT kits 
and found that, in our hands. Superscript works 
best for low- abundance non-coding RNA 
detection. 

* Note 2: We recommend using the same amount 
of RNA for each RT reaction. We typically use 1 
(ig total RNA for a reaction in 10 (iL. 

* Note 3: No-RT control reactions can be performed 
by making RT mixtures as describe above but 
without reverse transcriptase. This will ensure 

no false qPCR signals from DNA contamination 
in RNA samples. 

4. Perform qPCR by using SYBR® Green Master Mix. 
PCR for scaffold RNA and ACTIN2 should be 
separated. 

Step I: 95°C 3 min 
Step II: 95°C 10 sec 
Step III: Tm 30 sec 

Step IV: go to Step II 39 times + plate read 
(totally 40 cycles) 

Step V: increase from 65°C to 95°C, increment 

0. 5°C for 5 sec, plus plate read 
Step VI: End 

* Note 1: Because cDNAs are synthesized by using 
gene-/scaffold RNA-specific primers, PCR amplifi- 
cation that follows has to be nested PCR, 

1. e., positions of the qPCR primers on the 
DNA coding strand should be upstream of the 
gene-/scaffold RNA-specific RT primers. 

Note 2: A single peak in the melt curve 
generated by Step V should indicate the specificity 
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of the primers. This can also be examined by DNA 
gel electrophoresis. 
5. Use ACTIN2 as internal control to analyze the 
results (Refer to Steps V-11 to V- 17 for details). 

* Note: By starting with equal amount of RNA, 
different samples should show similar Ct value for 
ACTIN2 in qPCR reactions. Alternatively, another 
housekeeping gene can be examined to ensure 
similar RNA quality. 

<lll-2> Quantitative measurement of scaffold RNAs with 
nuclei fractionation 

1. Prepare Honda buffer (0.4 M Sucrose, 2.5% Ficoll, 
5% Dextran T40, 25 mM Tris-HCl, pH 7.4, 10 mM 
MgCl2, 0.5% Triton X-100, 0.5 mM PMSF, 0.1% 
pepstatin A, 0.1% aprotinin, 10 mM |3- 
mercaptoethanol, 8 unit/mL RNaseOUT) 
(modified from [9]). 

* Note: add PMSF, pepstatin A, aprotinin, 
|3-mercaptoethanol, and RNaseOUT before use. 

2. Grind 3 gram ~12-day-old plants with liquid 
nitrogen. 

* Note: We recommend using the same amount of 
plant materials for each sample. 

3. Transfer the power to a 50 mL Falcon tube prefilled 
with 15 mL cold Honda buffer. Homogenize by 
vortex vigorously. Keep the tubes in ice for 15-30 min. 

4. Filter the homogenate through 4 layers of Miracloth, 
followed by centrifugation at 1500 g at 4°C for 5 min. 

5. Re-suspend the pellet in 1 mL Honda buffer, then 
transfer to a new 1.5 mL eppendorf tube, followed 
by centrifugation at 1800 g at 4°C for 5 min to pellet 
the nuclei. 

* Note: Save the supernatant for cytosol analysis if 
wished. 

6. Gently rinse the nuclei pellet with ice-cold wash 
buffer (IX PBS, 1 mM EDTA, 8 unit/mL RNaseOUT). 

7. Re-suspend the pellet in 0.5 mL prechilled glycerol 
buffer (20 mM Tris-HCl, pH 7.9, 75 mM NaCl, 0.5 
mM EDTA, 0.85 mM DTT, 50% glycerol, 0.125 mM 
PMSF, 0.1% pepstatin A, 0.1% aprotinin, 10 mM 
|3-mercaptoethanol, 160 unit/mL RNaseOUT) by 
gently flicking the tube. 

8. Add 0.5 mL cold nuclei lysis buffer (10 mM HEPES, 
pH 7.6, ImM DTT, 7.5 mM MgCl2, 0.2 mM EDTA, 
0.3 M NaCl, 1 M Urea, 1% NP-40, 0.5 mM PMSF, 
Proteinase inhibitor cocktail, 0.1% pepstatin A, 0.1% 
aprotinin, 10 mM |3-mercaptoethanol, 160 unit/mL 
RNaseOUT). 

9. Gently vortex for 2 x 2 sec, incubate on ice for 
2 min, centrifuge at 14,000 rpm at 4°C for 2 min 
to separate the chromatin fraction and the 
nucleoplasmic fraction. 



10. Treat the supernatant (nucleoplasmic fraction) with 
proteinase K at 37°C for 1 hr, followed by traditional 
acidic phenol/chloroform extraction and ethanol 
precipitation for RNA. 

11. Gently rinse the chromatin pellet twice with 0.5 mL 
cold wash buffer and then dissolve in 1 mL Trizol 
for RNA extraction. 

12. Refer to Section III-l for RT-qPCR of scaffold RNAs. 

Additional files 



Additional file 1: Figure SI. Examples of methylation-sensitive Chop 
PGR and Chop qPCR. DNA methylation at a group of transposon loci 
were examined by Chop PCR or Chop qPCR. Mutants were compared 
with the wild type Arabidopsis [Col-O). For BsmF I, BstBI, or Hae III digestion, 
a DNA sequence from SKPl without the restriction sites was examined as 
loading control. For HpyCH4 IV, 5TE27040 was amplified as non-digestion 
control. [A] Chop PCR. Restriction enzymes BsmF I, HpyCH4 IV, BstB I, and 
Hae III are indicted on the right. The examined loci are annotated on the 
left. [B] Chop qPCR. Hae III was used for digestion. Methylation levels, as 
indicated by qPCR signals, in the mutants were relative to those in Col-O. 
SKPl was used as loading control. Error bars indicate SD, n > 3. Primers are 
listed in Table S3. Figure 52. Detection of Arabidopsis miRNAs by using 
TaqMan Small RNA Assay. Aerial portions of plants grown with (P deficient) 
or without (control) phosphate deficiency stress were examined. Each RT 
reaction used 80 ng total RNA and snoRWl was used as internal control. 
Error bars indicate SD, n = 3. Figure S3. Quality control of nuclei fractionation. 
qPCR detection of genomic DNA (gDNA) in the chromatin-associated fraction 
(pellet) and the chromatin-free fraction (supernatant) after nuclei fractionation. 
SKPl gDNA levels in the chromatin-free fraction were presented relative to 
those in chromatin-associated fraction. Error bars indicate SD, n = 3. Primers 
are listed in Table S3. (PDF 363 kb) 

Additional file 2: Table SI. Methylation-sensitive restriction enzymes for 
Chop PCR. Table S2. Arabidopsis 24nt siRNAs detected by using TaqMan 
Small RNA Assays. Table S3. Primers used for PCR/qPCR in this report. 
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